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Relative to the huge body of theory of linear time-invariant systems, very 
little of a general and precise nature is known about the network-theoretic 
properties of transistor circuits operating under large-signal conditions. 
One basic property P which a transitor network might have is that if the 
input approaches a constant, then the output approaches a constant which is 
independent of the initial conditions. In this paper we prove a stability 
theorem concerning a nonlinear differential equation that governs the 
behavior of a large class of networks. A corollary of this theorem asserts 
that if a certain condition is satisfied, then property P holds. 

We consider also the problem of estimating the rate of decay of transients 
in transistor networks and we prove theorems which allow us to make some 
often quite conservative, but definite, statements concerning limitations on 
switching speeds. A practical example considered shows that in some cases 
the bounds, which are frequently very easy to evaluate, can be quite useful. 

The proofs depend in an interesting way on the relationship between 
the static diode characteristic and the nonlinear capacitance associated with 
a semiconductor junction. 

I. INTRODUCTION AND DERIVATION OF THE DIFFERENTIAL EQUATION 

We initially consider the network of Fig. 1, which contains transis- 
tors, linear resistors, voltage sources, and current sources. Each 
transistor is represented by a model of the type shown in Fig. 2 (see 
Gummel 1 and Koehler 2 ) which takes into account nonlinear dc proper- 
ties as well as the presence of nonlinear junction capacitances. Asso- 
ciated with this model are six parameters: «/, «,-, t c , t c , c c , and c c (all 
positive constants; «, < 1, «, < 1^» and two nonlinear functions /<>(•) 

and / c (-)- 

Concerning / e (-) and f {-), for our purposes it is necessary to as- 
sume only that 
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TRANSISTOR I TRANSISTOR 2 



Fig. 1 — General network containing transistors, sources, and resistors. 

Assumption 1: For each transistor: /<,(■) and f c (-) are strictly -mono- 
tone increasing mappings of the real interval (—00, 00) into itself; 
/ e (0) = / c (0) = 0, and / e (-) and / c (-) are continuously differentiable 
on ( — 00, 00). 

The functions /<,(•) and /„(•) of Gummel's model 1 are of simple 
exponential type and satisfy Assumption 1. 
From Fig. 2: 

i. = J t kft + TJM] + /.(!>.) - «r/.(0. 
t. = f t [C C V C + TJM] ~ a,1.(v.) + M. 

Suppose that the network of Fig. 1 contains p transistors; for k = 
1, 2, ' " , p, let ^2fc-i and v 2 k, respectively, denote the emitter to base 
voltage and the collector to base voltage of the fcth transistor. Simi- 



a r t. 



atf Of 




Fig. 2 — Transistor model. 
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larly, for k = 1, 2, • • •, p, let u k -i and u k , respectively, denote the 
emitter current and the collector current of the fcth transistor (with 
reference polarities as indicated in Fig. 2). Then, with v = (v lf v 2 , 
• • • , v 2p ) "', i = (h ,12, • • • , i2p) "\ 1-ik-\ ( • ) and c^-\ the / e ( • ) and c of 
the fcth transistor, and / 2fc (-) and c 2fc the /„(•) and c c of the fcth tran- 
sistor, 

i = J ( [C(v)] + TFW U) 
where, for ;' = 1, 2, . . . , 2p, 

[C(p)l = C iVj + r,/,(l>,) (2) 

[F(v)l = ffa), (3) 

and T = T, © T 2 © • • • © !T P , the direct sum of p 2 X 2 matrices 
T t in which 



T k = 



1 -«; 



(A) 



(*) 



for fc = 1, 2, • • • , p. 

We assume that the linear resistive portion of the structure of Fig. 1 
introduces the constraint 

i = -Gv + B (4) 

in which G is a conductance matrix and B is an element of the set (B 
of all real bounded continuous 2p-vector- valued functions of t on [0, <»). 
From equations (1) and (4) 

j t [C(v)] + TF(v) + Gv = B. (.5) 

Let u = C(v). Since all of the c,- and t, are positive, and each of the 
/,-(■) is continuous and monotone increasing, there exists aC" (•) such 
thatv = C~\u). Thus, 

j t + TF[C-\u)] + GC'-\u) = B. (6) 

The Jacobian matrix J u of TF[C'\u)] + GC~\u) is 

T diag { JjHf] J + G diag ( M * . J 
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in which for all j = 1,2, • • • , 2p 

with each of the </,•(•) continuously cliff erentiable. 

Since J u is continuously dependent on u, and || J u || (|| • || any norm) 
is bounded from above uniformly in u, it follows that there exists a 
constant L such that 

|| TF[C~\u a )] + GC~\u a ) - TF[C- l (u b )\ - GC~ l {u b ) \\ 

^ L || u a — u b || (7) 

for all u a and u b belonging to real Euclidean 2p-space E 2 ". In particular, 
we have 

|| TF[C-\u)] + GC-\u) - B || £ L || u || + ||£ || (S) 

for all t ^ and all u e £ l2p . Therefore (see, for example, Nemytskii and 
Stepanov 3 ), for any initial condition u £ E 2p , there exists a unique con- 
tinuous 2p-vector-valued function u( • ) such that w(0) = v and equation 
(6) is satisfied for all t > 0. In other words, under the assumptions we 
have introduced, it makes sense to study the properties of the solution 
of the equation 

j t + TF[C-\u)\ + G[C-\u)} = B, / ^ K0) = u ] (9) 

II. STATEMENT OF RESULTS, AND EXAMPLES 

We need the following definitions. 

Definition 1: A real matrix M of arbitrary order n is strongly col- 
umn-sum dominant if and only if for all j = 1, 2, . . . , n 

•»« — S I m a I > 0- 

An important property of T is that it is strongly column-sum dominant. 
Definition 2: We shall say that a real matrix M of order 2p is an element 
of 2D if and only if there exists a diagonal matrix diag (d, , d 2 , • • • , d 2v ) 
with each d f > such that 



(A) ^ "2A-1 ^ 1 
rf 2 A «r 



a., < —z < - (n 



for A; = 1,2, • • ■ ,7?, and diag {dx , d 2 , • • ■ , c/ 2p ) Af is strongly column- 
sum dominant. 

Our main result* concerning equation (9) is: 

* Proofs of all results iu this section are given in Section III. 
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Theorem 1: If G t 2D, and u a (-) and u b (-) satisfy 
*■ + TF[C'-\u a )] + G[C-\u a )) = B a (t), 



t ^ 



(10) 



jf + TF[C~\u,)] + G[C'-\u b )] = B b (t), 



t > 



(ID 



with B a z®> and B b e 03, a??cZ i/ [B a (t) - B b (t)] -» (Me zero wctor o/ S 2 ") 
ast—> oo , Men [w a (0 — u b (t)] -+ 9 as t ^> » . 
An interesting corollary of Theorem 1 is 

Corollary 1: Referring to equation (9), if G e 3D, and ?7 Mere ernsis 
a constant vector B*, such that [B(t) — B K ] — » as i — > oo, Me?i Mere 
exists a constant vector u x such that [u(t) — u K ] —> 9 as t —> oo , a/?rf m m t's 
independent of the initial condition w„ • J" particidar, if B v = 9, then 
Moo - fl- 
it is interesting to observe that Ge3) whenever the base leads of all 
transistors are connected together and there is a resistor between the 
emitter and base, and between the collector and base, of every transistor, 
for then G is strongly column-sum dominant. Also it is easy to give 
examples of conductance matrices which are not strongly column-sum 
dominant, and which belong to SD. For instance, for the network of 
Fig. 3. 




A/VV 



♦t 



1 
J 



9a >0 
9 b >0 



Fig. 3 — Single-transistor network. 



G = 



Qa + 06 ~g t 

— Qb Qb 



and diag (o\ , d 2 )G is strongly column-sum dominant for d 2 
some d x such that 

1 * 
a t < d t < — 



= 1 and 



* More generally, G of order 2p with positive diagonal elements belongs to D 
whenever it is possible to obtain a strongly column-sum dominant matrix from 
G by adding an arbitrarily small positive quantity to a single diagonal element. 
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Illl 



to 



A/W- 
IA 




IA 



X. 




Fig. 4 — A two-transistor circuit. 

As another example, consider the circuit of Fig. 4, for which 

473 -10 10 -11 

-10 473 -11 10 

10 -11 11 -10 

-11 10 -10 11 



G = n 



Since diag (1, 1, 22, 22)G is strongly column-sum dominant, G t ID. 
Finally, for the network shown in Fig. 5, 



G_ 21 



11 -10 10 -11 

-10 11 -11 10 

10 -11 11 -10 

-11 10 -10 11 



In this case, G is obviously singular and hence does not belong to 2D. 
Suppose that the source current of Fig. 5 i (t) is a constant and that the 
transistor functions /i(-), U( m ), / 3 (-)> and /<(•) are all bounded from 
below by the constant b (this is certainly an assumption consistent with 
our earlier assumptions and with the character of transistor models 

in 

m 




AA/V 
in 




Fig. 5 — Transistor circuit for which the dc equations may have no solution. 
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ordinarily used.) We wish to show that here for sufficiently small i , 
there does not exist a constant vector u„ such that [u(t) — w»] — > 6 
as t — > co . 

Suppose that u(t) -»«,,a constant vector, as i — » °° . Then there 
would exist a 2p-vector y M such that u„ = C(v„) and 

TF( Va3 ) + Gi;. = B 

with £ = (z„ , 0, 0, • • • , 0) tr . Let rj denote the 2p-row-vector (1, 1, 1 , • • • , 1) . 
Then 

v TF(v M ) + r,Gv a = tjB. 

But r)Gv M = 0, and hence 

which does not possess a solution «;„ if 

i < b i; [i - «n + [i - aj»]. 
*— i 

2.1 Estimation of the Rate of Decay of Transients 

Theorem 2: If the hypotheses of Corollary 1 are satisfied with B(t) = B^ 
for t ^ 0, then 

E <*, I «i(0 - «-/ I ^ exp (-£*) E d, 1 1*,(0) - it., |, t* 



)-i 



i-i 



/or euery se£ o/ positive constants d x , d 2 , • • • , d 2v such that 

< K = min min {— (1 - ^d? 1 **,), - (ff„ - E <W I 0« \)\ 

i \Ti c i i*i > 

in which — a { is the nonzero off-diagonal term in the jth column of T, 
and 3, = d i+1 for j odd and U t = d,_i for j even. 

It is easy to show that G c £) imphes that there are positive constants 
di , j = 1, 2, • ■ • , 2p, such that %. > 0. 

As an example of the application of Theorem 2, consider the problem 
of estimating the switching time of the single-transistor inverter circuit 
of Fig. 6 in which a, = 0.968, c. = 2 X 10" 12 fd, t, = 1.7 X 10" 10 second, 
a r = 0.583, c e = 1.7 X 10" 12 fd, and r c = 2.62 X 10" 8 second. Here 
(in mhos) 

1.1886 X 10" 3 -1.01215 X 10" 3 



G = 



-1.01215 X 10" 3 1.01215 X 10- 
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+ 
V2 



IOK 

I — wv 



«o 




-=- 10V 




Fig. 6 — Practical logical-inverter circuit. 



which takes into account a buDc base resistance of 280 ohms and a bulk 
collector resistance of 18 ohms. The circuit is initially at steady state 
with e(t) = 0.3 volt for t < 0. For t ^ 0, e(t) = 10 volts, and as t -> °o, 
u(t) — > We , some constant vector. With d 2 = 1, the number K is the 
smallest of the four quantities: 0.58(1 - 0.968O X 10 10 , 0.5(1.1886 - 
1.012150 X 10 9 , 0.3815(1 - 0.583d,) X 10 8 , and 0.58(1.01215) 
(1 - ch) X 10°. 

It is clear that d x must satisfy 0.968 < d x < 1 in order that K > 0. 
Then optimal choice of d r (that is, the choice that yields the largest 
value of K) is approximately 0.9709. For d x = 0.9709, K = 1.66 X 10 7 . 
Let the "charge switching time" t, denote the smallest value of t such 
that ^y„i | u t (t) — Wooi | is less than or equal to two percent of y^„, 
| w,(0) — Wo,,- | for all t ^ t, . Then our upper bound on t, is approxi- 
mately 4 X (1.66)" 1 X 10~ 7 w 241 nanoseconds. The actual value of t, , 
as determined by numerically integrating the system of two nonlinear 
differential equations is approximately 57 nanoseconds. Thus, for this 
circuit, Theorem 2 provides a very easily evaluated and useful upper 
bound on t. . 

Finally, we state a result which provides an often rather conserva- 
tive but easily evaluated lower bound on the rate of decay of tran- 
sients. 

Theorem 3: With B a constant real 2p-vector, let 

j t + TF[C-\u)\ + GC-\u) =B, t£0. 

If there exists a constant 2-p-vector u„ such that [u(t) — u„] — > 8 as t — > °o , 
then for any choice of positive constants dj , j = 1,2, • • • , 2p: 



2p 

z 

f-I 



2 d, | w,(f) - Woo,- | ^ exp (-Kt) 2 d f | w,(0) - u xj |, t ^ 
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in which 

R = max max 



{-(l + flrfT 1 *,),*- jtdtf | fir., [V 



wAere — a,- is the nonzero off-diagonal element in the jth column of T, and 
3j = d,+i for j odd, and Sj = d f _j for j even. 

The arguments used to prove the results stated in this section can 
be modified in a straightforward manner to prove far more general 
results concerning networks that contain diodes, capacitors, and in- 
ductors, in addition to the elements of the structure of Fig. 1. Some 
of these more general results are described in Section IV. 

III. PROOFS 

3.1 Proof of Theorem 1 
We first show that 

F[C-\u a )] - F[C-\u b )] = D 1 (u a - u b ), t ^ (12) 

and 

C~\u a ) - C-\u b ) = D 2 (u a - u b ), t ^ (13) 

with Di and D-, diagonal matrices dependent on t and possessing some 
special properties. 

For ; = 1, 2, ••• , 2p, let Qj{u a i) = [C" 1 (««,)]> and gj(u bj ) = 
[C _1 (ii h ) ],-. Then, using equation (2), 

u ai - u bj = Cy[0,.(w o/ ) - Oi(u bi )] + T,{/ J -[ff,(w a ,-)] - fAgi(u bi )]}. 
Thus if u^ ** u bi j 

1i[9i(Uai)] - f l [g i (u b ,)] = r,(u ai , u bi ) 

U ai - U bi C,- + TiTffai , u bi ) ' 

in which (for u ai j* u bi ) 

r (,. ,. N _ Ji MUaj)] ~ 1iMu b ,)] 

r A u aj , U bi ) — . r - , . 

9i\U ai ) - QAU bi ) 

In a similar manner we find that for all u ai 9^ u bi : 

g/ftO ~ Qifabi) = 1 

u ai - u bi Cj + Tir(u ai , u bi ) 

Now, let us define for j = 1, 2, • • • , 2p 

ri(u ai , u bi ) = f'i[gi(u ai )] 
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when u ai = u bi . Then since u ai and u bi are continuous on [0, oo), it 
follows (see Appendix A) that r t (u ai , u bi ) is continuous on [0, ■»). 
Since rj(u ai , u bj ) is nonnegative, it is clear that both 

rjfaai , Ubi) 



c,- + rpfyln , u bi ) 
and 

1 

C,- + T,T,-(M a ,- , U bj ) 

are continuous on [0, oo). Moreover equations (12) and (13) are satisfied 
with 

A = diag / r,{u ai u b ,) \ 

& \c,- + T^fai , u bj )) ' 

D 2 = diag \—r \ r}« (15) 

At this point we have 

j t (u a - u b ) + (TD t + GD 2 )(u a -u b ) = B a - B b , t^O (16) 

with TDi + GD 2 continuous on [0, °o). 
We need the following lemma. 

Lemma 1*: Let M(') be a continuous real n X n matrix-valued junc- 
tion of t defined on [0, <») such that there exist positive constants e and 
Ci , c 2 , • ■ ■ , c„ , with the property that for j = 1, 2, • ■ • ,n and allt^O 

™>a - £ Wj l | m {i | ^ e. 
Let xbe a differentiate real n-vector-valued function on [0, °o ) such that 

^ + Mx = 0, * 2> 0. 

dt 

Then there exists a constant k such that for i = 1, 2, • • • , n, and allt ^ 

\xS) | g k exp(-rf). 
Moreover, k depends only on the c, and //ie initial values .r,(0). 



* In Ref. 4, Rosenbrock states a similar result, but does not give a rigorous 
proof. He considers the case in which c$ = 1 for ; = 1, 2, • • • , n. 
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Proof of Lemma 1: Let the functional s be denned in terms of an 
arbitrary continuously differentiable scalar function <?(•) by 

8(<p)(t) = 1 if <p(t) > or if <p{t) = and tp'{t) > 

= -1 if tp(t) < or if <p(t) = and tp'ify < 

= if tp{t) = and <p'(f) = 0. 

Then for t ^ 0, 

E cMxMxW - - E c,s(.r,)(0 E »**«« 

i i 1 

= - E z y E c,s(a;,)(0m, y 

= - E xfiifoMitu - E x > E e<a(*O(0i»« 

£ - 2 c,wi,- ; | as, | + Si ** I E <>i I m„ I 

^ -«El<***|. 
/ 

7 + 

But E- c*«(**)(0*< is ec l ual to j? E» I c i x i li tne right-hand derivative 

of E» I c * x * I t see Appendix B; the derivative of | x f | need not exist 
at points i at which Xj(t) = 0]. Therefore 

^Ek*,l ^ -ell^l, t^O 

from which it follows that 

E I <*r,(*) | =s exp (- c E I c,x,(0) I, I £ 0. D 

i i 

If Af(-) satisfies the conditions of Lemma 1, then it is easy to show 
that the unique continuously differentiable n X n matrix-valued function 
X denned on [0, » ) which satisfies 

—■ + MX = 0, / £ [X(0) = /] 
at 

possesses the property that (for any norm || • || on E") there exists a 
constant K x such that 

|| X(t)X(r)- 1 1| ^ K, exp [-«(«- r)] 

for all « ^ r. 

Returning now to equation (16), assume that [TD t + GZ) 2 ] satisfies 
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the conditions on M(-) of Lemma 1. Then with Y the solution of 
|y-+ [TD t + GD 2 ]Y = 0, t^Q [Y(0) = I] 

we have 
u a (t) - u b (t) = Y(t)[u a (0) - u b (0)] 

+ f Y{QY(t)- 1 [B.(t) - B b (r)] dr, t ^ 0. 

Jo 
Therefore, for t ^ 

|| u a (tj - ih(t) || ^ || Y(t)[u a (Q) - uM]\\ 

+ f || Y^Yir)- 1 \\-\\B a (r) -B b {r) \\ dr 
Jo 

^ || F(/)K(0) -u b (0)] II 

+ £ a f exp [-«(! - r)] || B.(r) - 5 6 (t) || dr 

for some positive constant K 2 . Since || B a (r) — B h {T) \\ — > as t — > oo , 
it follows that || w„(i) — i*(,(f) || — » as £ — > «> . 

It remains only to prove that [TZ)i + GZ) 2 ] meets the conditions 
imposed on M(-) of Lemma 1. Since Gt3), there exists a diagonal 
matrix diag {cU , d 2 , ■ • • , d 2v ) with d f - > f or j = 1, 2, • • • , 2p and 



'/ < — -, — < ~u 

«2fc «r 



for /c = 1,2, • • • , p such that both 

diag (di ,d 2 , ■■• , d 2p )G 
and 

diag (d 1} d 2 , ■■■ , d 2p ) T 
are strongly column-sum dominant. Thus for j = 1, 2, • • • , 2p 

'„ - Z «MT' I '</ I > o 
!7,/ - ZddJ 1 \9u I >o. 
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Let W - TD t + GD 2 . Then, for j = 1,2, ■■■ , 2p, 

r 1 

w >> = l a „ ■ ' + On 



c, + r-r, »" c s + T i r i 
and 



T. d,dj l I w« I = S WT 1 



iVi 



'■•/ 1 r : - 1 - 9a 



c, + r,r, y " c, + T,T,- 



Therefore 



w« - E «W I "« I ^ r^r— &, - £ wr 1 1 t« I) 

■Vj W I TjTj ,V; 

Since r,- ^ 0, the right side of equation (17) is bounded from below by 
some positive constant e uniformly in t and j. □ 

3.2 Proo/ o/ Corollary 1 

By Corollary 3 of Ref . 5 there exists a unique v t E 2p such that 

TF{v) + Gv = B„ (18) 

whenever G is such that all principal minors of T~ G are positive. 
In Reference 5 it is proved that T~ l G will have this property if T~ l G can 
be written as A~\B with both A and B stongly column-sum dominant. 
Let H = diag (d x , d 2 , ■ ■ • , d 2p )G be strongly column-sum dominant 
with all dj > and 



U) ^ d-2k-l ^ 1 
«2fc «r 






for A; = 1, 2, • • • , p. Then £/ = diag (d, , d 2 , ■■■ , d 2p )T is strongly 
column-sum dominant, and T~ l G = U~*H, which proves that equation 
(18) possesses a unique solution v. 

With v the solution of equation (18), let w w = C(v). Clearly if B a = 6, 
then Woo = 0- Let u b satisfy 

^ + TF[C-\u h )] + G[(T\u h )] =B m , t^O 

with u b (0) = Ma, . Of course, u b (t) = u x for all t ^ 0. By Theorem 1, 
[w(£) — Moo] — » as £ — » o° , independent of w . 
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3.3 Proof of Theorem 2 

Following the proofs of Theorem 1 and Corollary 1, 

4-Au- «J + (m + GD 2 )(u - O = 0, / ^ 
at 

in which 

r,(Uj , Woo,) 



and 

1 



D, = diag 
D 2 = diag , 

a {C, + T,T,(Wy , M M) ) 

Therefore 

%T,di\ u,{t) - «„,- | £-££«*# I »/(0) - «-< I . * fc 
in which 
£ = min min S— (£,-,• - £ d.d" 1 | «,-, |), - (g if - £ <W I 9u l)r 
But for j = 1, 2, ••• ,2p 

«,-i - E dwr 1 1 tu 1 = 1- Wr'«/ • □ 

3.4 Proo/ o/ Theorem 8 

Since Ti^C" 1 ^)] + GC~\z) depends continuously on z t E 2p , u„ 

satisfies (see Ref. 6) 

TF[C~\uJ] + GC-\u„) = B. 

Therefore, following the proofs of Theorem 1 and Corollary 1, 

j t (u - u.) + (TD, + GD 2 )(u - O = 0. /• = 



in which 

V,- - 

D 2 = diag 



„, , uia ^ r,(t, f H,,) 1 

1 7^^ , U xi )) 



and 

1 
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For any z c E 2v , let || z || denote £ d f \ z i I • Then > for l = ° 
j t (u-u a ) =|| (TZ), + GZ) a )(i* - w.) II 

= m ? X \ (1 + ** a '° c,. + r,r,(«, , «.,) 

+ t ddl 1 I 0<, I - j—kr. ~ } I MO - «- II- 



But, since r y (w,- , w w ,) ^ 0, 

2p 1 



Thus 



5 («-m0 



< 7? 1 1 M - M. 



/ > 0. 



Clearly, 
d , N 

5 («-«-) 



= lim - II u(t + e) - w- - w(0 + w 



(19) 



/ 2: 



Also, for £ ^ 0, the limit 

lim - [|| u{t) - u x || - || u(t + e) - w M ||] 

exists and is equal to - — || u - u. || in which as before — denotes 

the right-hand derivative (see Appendix B). But, since for any e > 
and t = 0, 

|| u(t) - u, || - || u(t + e) - u x || ^ || u(t + e) - u„ - u(t) + «. ||, 
we have 



(ft 



U — U a — : 



It (" - «") 



Therefore, using equations (19) and (20), 

^|| tl - ^ || = -it \\u-u 



t > 0. 



/ = 



(20) 
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and, for t ^ 0, 

II u - W. || ^ exp {-lit) || w(0) - u x ||. D 

IV. A SIGNIFICANT EXTENSION 

We can easily extend our results to cover an interesting class of 
networks containing diodes, capacitors (not necessarily linear) , and 
(not necessarily linear) inductors, in addition to the elements of the 
Fig. 1 network. 

Let each diode be represented by a model of the type shown in 
Fig. 7 in which 

id - -g [c d v d + T d j d (v d )} + f d (v d ), 

with c d and r d positive constants. Assume that /,,(•) satisfies the con- 
ditions placed on /.(•) and / e (-) of the transistor model. Let there be 
q diodes and let v 2p+k and 4»+* (k = 1, 2, • • • , q) be the voltage and cur- 
rent associated with the fcth diode. 

Suppose that the /cth capacitor (we assume that there are r capacitors) 
is governed by 

"77 [C2p + g + k\P2j> + a + k)i = ^2p+a + A; 

for k = 1, 2, • • • , r, where c 2j>+Q+k (-) is a strictly-monotone-increasing 
continuously-differentiable mapping of E 1 onto itself such that c 2 p +a+4 (0) 
= and the slope of c 2p+ „+*(•) is uniformly bounded from above and 
from below by positive constants. 

Finally, let there be s inductors which introduce constraints 

~T. l'2j> + o + r + H^2ji + fl + r + i/J = ^2p + o + r + t 




Fig. 7 — Diode model. 
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for k = 1, 2, • • • , s, in which each Z 2p+fl+r+i (-) is a function of the same 

type as the c 2p+ „+*(•)• 

Assume that the linear resistive portion of the network introduces 

the constraint 

i = -Hv + B, Be(B 

in which i = (t\ , i 2 , • " , *a+»+«+n i'2p+<,-rr+i , ■ ■ ' » v 2v + Q+T +,) , v = 
(vi , v 2 , ■■■ , y 2p+a+r , i 2 p +fl+ r + i , • • • , *2p + , + r + .)", and H is a constant 
hybrid-parameter matrix of order (2p + g + r + s). Then 

| [ff(f»] + 7'F(f) + Hv = B 

where 

[£(»)], = [CW1#, *= 1,2, ••• ,2p 

= Cfi, + r,/,fr#), i = °-V + 1, 2p + 2, • • • , 2 V + g 

= «#(»#), ; - 2p + g + 1, • • • , 2p + g + r 

= *,(«#). i = 2p+g+r+l, ■•• ,2p + g+r + s; 

f is the direct sum of matrices T © I a r+ . , in which J, is the identity 
matrix of order q and r+ . is the zero matrix of order (?• + s), and 

[F(v)l = [TOL , J = 1,2, ••• ,2p 

- />>), j = 2p + l, ... ,2p+g. 

Under our assumptions C(-) -1 exists and, with & = (7(0), 

^ + fF[C-'(u)] + fffl- l (fl) - -B. (21) 

Let 2& denote the set of all real matrices M of order (2p + g + r + s) such 
that there exist positive constants d Y , d 2 , • ■ ■ , d 2j>+q+r+ , with the 
property that 

U) ^ (hk-i ^ 1 
ctf < ~~~ ; <. <*) 
a 2k a T 

for fc = 1, 2, • • • , p (when p^O) and diag (d, , d 2 , • • • , d 2ptatr+ .)M 
is strongly column-sum dominant. 

With straightforward modifications of the arguments already pre- 
sented, we can prove (i) that for each u Q t _B (2p+0+r+ * ) equation (21) 
possesses a unique solution defined on [0, °°) such that w(0) = u , and 
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(it) the analogs of Theorems 1, 2, and 3 and Corollary 1. To be more 
specific, the analogs of Theorem 1, Corollary 1, and Theorem 2 are: 

Theorem 1' : If H e £>, and u a and u b are solutions of equation (21) 
with B = B a and B = B b , respectively, for t ^ 0, and if [B a (t) — B b (t)] 
— » 6 [the zero vector of 2£ <a,+ « +r+0 ] as t — > « wz*A £ e 63 and 5 6 e (B, then 
[u a {t) — u b (t)] —> 6 as t —> oo. 

Corollary 1' : Referring to equation (21), if H z SD, and i/ f/iere ezzs£s 
a constant vector B„ such that [B(t) — B x ] —> 6 as t — > oo , then there exists 
a constant vector u n such that [u(t) — u K ] — > 6 as t — > °o , and u M is in- 
dependent of the initial condition u . In particular, if B„ = 6, then u„ — 6. 

Theorem 2': If the hypotheses of Corollary 1' are satisfied with B{t) = 
Boo for t ^ 0, then with j = (2p + q + r + s), we have 

£ d, | Ui(t) - u ai | ^ exp {-Kt) £ d t | uM - u ni |, t £ 

i-i y-i 

/or ever;/ set of positive constants d, , d 2 , • • • , d 2p+Q+r +, such that < K = 
min {Ki , K 2 , K 3 ] where 

K x = min min \— (1 - cljdj 1 *,), - (fl,-,- - £ d.d^ 1 I On \)\ 
K 2 = min min \— , — (g ti - £ d t dj l \ g if \)\ 

2p+lS/S2p+« k?"j W »Hj J 

K 3 = min {f ( 9ii - £ d.dr 1 I 9u I)} 

in which s,- = sup c'j(-) for j = 2p -\- q -\- I, ■ • ■ , 2p -\- q -\- r; s,- = 
sup l'j(-) for j = 2p + g + r + 1, • • • , 2p + o + r + s; — a,- is Me 
nonzero off-diagonal term in the jth column of T; and a 1 / = d i+x for j odd 
and d>i = d,-_! for j even. Moreover there exists one such set of constants {d f ). 

V. FINAL COMMENTS 

The results presented here are quite encouraging in that they are 
concerned with the equations of reasonably realistic nonlinear network 
models, and provide some understanding of a precise nature in an 
area where there is a great need for many results of similar type. 
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APPENDIX A 

Proof that r,(w„,- , u bi ) is continuous. 

It is clear that r,(u , , u bi ) is continuous at each point t such that 
Uaj (t) 9± u bi (t). Suppose now that t is such that u ai {t) = u bi (t), and 
let e > be given. Since u ai , u bi , g and /; are continuous, there exists 
5, > such that 

\j'i MucS + v)]) -ft Mu*M}\ ^ * 

for all | v I ^ &i • Then for | i?| = fi i either u ai {t + r,) = u bi (t + r?) in 
which case 

I r,[u a1 (t + t,), u bi {t + tj)] - r,[« ai (Q, u bi {t)\ | ^ <= 
or u bi (t + j?) ?* u bi (t + v) and (using the mean-value theorem) 

r,[u.,(« + *), «,,« + u)] - ^^ + ,)] _ ^(j + ,)] 

= rm 

in which 

I £ - ft[ti»f(0] I ^ max { | g,[u a & + i?)] 

-fclMfll I , I fr[«w« + i)] - tf#K*(01 I }• 

In the latter case, there exists 8 2 > such that | /'(£) - f lff/[«-*(0]} I 
^ e for all | rj | ^ 5 2 . Thus for all | i? \ ^ min { 5, , 5 2 } , we have 

I r,[u a ,(t + v), «m(* + *)] - »K*(0i «w(*)] I ^ 6 - 

APPENDIX B 

Proo/ too* the Right-Hand Derivative oj | x t | exists and is equal to 
If « is a point such that 3,(0 ^ 0, then it is clear that 

At t such that x,(t) = ° and xj(0 ^ 0, 



54 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 19C9 

Finally if x f (t) = and x' § (t) = 0, then 

= lim | xm | = lim I Xi(t + 6) I = ^ | x, |, 

t-*0+ «-»o+ e tt* 

tstit+t 

since x { is continuously differentiable. 
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